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1.  INTRODUCTION 


Quantum  chemistiy  is  an  invaluable  tool  in  determining  details  of  reaction  mechanisms  as  well  as  in 
characterizing  other  properties  of  molecules  (Schlegel  1987;  Hehre  et  al.  1986).  A  proper  description  of 
the  saddle  point  and  the  reaction  path  through  quantum  chemical  calculations  is  of  great  importance 
because  most  details  of  these  are  inaccessible  through  experiment  (Moeller  1934).  Theory  can  give 
stractures,  characteristics,  and  energies  of  species  along  the  entire  reaction  path  through  transition  state 
and  intrinsic  reaction  coordinate  (IRQ  calculations  (Schlegel  1987;  Hehre  et  al.  1986).  A  complete  set 
of  such  information,  however,  requires  computation  of  energies  and  energy  derivatives  of  many  points  on 
the  potential  energy  surface  (PES).  This  can  be  accomplished  at  the  Hartree  Fock  (HF)  level  for  many 
molecules  using  modest  computer  resources.  However,  HF  does  not  perform  well  in  predicting  reaction 
barriers  and  properties  of  structures  with  stretched  bonds.  Electron  correlation  must  be  introduced  in  order 
to  improve  the  accuracy  of  these  calculations.  One  of  the  least  expensive  of  the  ab  initio  electron 
correlation  methods  is  second-order  Moeller-Plesset  (MP2)  perturbation  theory  (Moeller  1934).  However, 
post-HF  treatments  such  as  these  scale  as  N^  or  greater,  where  N  is  the  number  of  basis  fimctions, 
effectively  limiting  the  number  of  calculations,  size  of  basis  set,  and  size  of  the  system  studied.  For 
systems  in  which  the  reaction  path  contains  a  multitude  of  transition  states  and  minima  and/or  a  large 
number  of  atoms,  a  complete  determination  of  the  reaction  path  becomes  prohibitive. 

Within  the  last  decade,  density  functional  theory  (DFT)  has  offered  a  computationally  less  expensive 
yet  reasonably  accurate  alternative  to  ab  initio  methods  for  including  correlation  corrections  in  calculating 
molecular  properties  such  as  geometries,  frequencies,  and  energies  (Hohenberg  and  Kohn  1964;  Kohn  and 
Sham  1965;  Zeigler  1991;  Labanowski  and  Andzelm  1991;  Politzer  and  Seminario  1995).  DFT  can  be 
applied  to  large  polyatomic  systems  since  the  computational  effort  of  density  functional  methods  scales 
as  N^.  Most  of  the  successes  of  DFT  have  been  in  descriptions  of  stable  molecules.  The  more 
sophisticated  nonlocal  DFT  (NDFT)  methods  provide  reaction  energies  and  molecular  properties  that  are 
typically  better  than  the  HF  method  and  are  closer  to  (and  in  some  cases  better  than)  the  correlated 
methods  (e.g.,  MP2)  (Porezag  and  Pederson  1995;  Broclawik  et  al.  1995;  Jursic  1995;  Deng  and  Zeigler 
1995;  Johnson  et  al.  1994;  Deng,  Branchadell,  and  Zeigler  1994;  Sosa  and  Lee  1993;  Fan  and  Zeigler 
1992;  Andzelm,  Sosa,,  and  Fades  1993;  Fan  and  zieigler  1990).  If  DFT  can  correctly  describe  all  of  the 
important  regions  of  the  potential  energy  surfaces,  such  as  transition  states,  stable  intermediates,  and  points 
along  reactions  paths,  detailed  reaction  mechanisms  can  be  extracted  for  systems  that  are  not  accessible 
either  through  experiment  or  ab  initio  studies.  Additionally,  due  to  its  modest  computational  requirements. 
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DFT  altars  to  be  the  most  promising  quantum  mechanical  method  to  incorporate  into  dynamics 
simulations  for  the  direct  calculation  of  forces.  If  DFT  can  be  shown  to  be  reliable  in  describing  chemical 
reactions,  one  of  the  main  limitations  of  dynamics  simulations  (the  use  of  analytic  functions  fitted  to  a 
limited  set  of  data  for  the  description  of  the  interatomic  interactions)  will  be  removed  for  many  cases. 
Therefore,  it  is  imperative  that  studies  such  as  those  presented  here  compare  the  performance  of  DFT 
predictions  of  transition  states  and  reaction  paths  with  well-documented  ab  initio  methods  such  as  MP2 
and  QCISD(T)  to  determine  both  advantages  and  limitations  of  the  methods. 

There  are  a  few  studies  that  have  apidied  DFT  to  transition-state  and  reaction-path  predictions  (Deng 
and  Zeigjer  1995;  Johnson  et  al.  1994;  Deng,  Branchadell,  and  Zeigler  1994;  Sosa  and  Lee  1993;  Fan  and 
Zeigler  1992;  Andzelm,  Sosa,  and  Eades  1993;  Fan  and  Zeigler  1990;  Baker  et  al.  1995).  The 
performances  of  local  and  NDFT  methods  when  compared  to  ab  initio  are  not  consistent.  In  some  cases, 
DFT  predictions  of  energy  barriers  are  as  good  as  ab  initio  or  in  closer  agreement  with  experiment  than 
ab  initio  (Deng  and  Zeigjer  1995).  In  other  cases,  DFT  fails  completely  in  correctly  describing  reactions, 
such  as  predicting  transition  states  as  stable  points  (Johnson  et  al.  1994),  or  as  ntmexistent  on  the  PES 
(Baker  et  al.  1995),  or  as  having  relative  energies  that  are  lower  than  the  minimum  to  which  it  cormects 
(Sosa  and  Lee  1993;  Fan  and  Zeigler  1992).  Also,  local  DFT  has  predicted  minima  on  a  PES  that  aj^ar 
to  be  artifacts  of  the  local  density  approximation,  and  do  not  exist  when  the  nonlocal  corrections  are 
incorporated  (Fan  and  Zeigler  1992).  In  general,  DFT  predicts  barrier  heights  that  are  lower  than 
experiment,  whereas  ab  initio  calculations  typically  overestimate  the  measured  values  (Porezag  and 
Pederson  1995;  Deng,  Branchadell,  and  Zeigler  1994).  The  inconsistencies  in  performance  of  DFT  for 
different  chemical  systems  indicate  the  need  for  metrics  such  as  those  presented  here  to  gauge  the 
usefiditess  and  suitability  of  DFT  methods. 

In  this  paper,  we  report  NDFT  calculations  for  two  possible  decomposition  mechanisms  of  sywi-triazine 

(H3C3N3): 


.^m-triazine - >  (HCN)3 - >  3HCN  (I) 

sym-triazir^ - >  HCN+  H2  C2  N2 - >  3HCN  (II) 

where  (HCN)3  denotes  a  hydrogen-bonded  cyclic  cluster  of  HCN  molecules,  and  H2C^N2  is  a  covalently 
bonded  dimer.  These  reactions  have  been  characterized  in  a  previous  ab  initio  study  (Pai,  Chabalowski, 
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and  Rice  1996).  In  Uiat  work,  critical  points  along  the  reaction  coordinates  for  (I)  and  (II)  were  located 
through  MP2  geometry  optimizations  using  the  6-31G**  (Hehre,  Ditchfield,  and  Pople  1972;  Hariharan 
and  Pople  1973;  Gordon  1980),  6-311-h-G**  (McLean  and  Chandler  1980;  Krishnan  et  al.  1980),  and 
cc-pVTZ  (Woon  and  Dunning  1993;  Kendall,  Dunning,  and  Harrison  1992;  Dunning  1989)  basis  sets  and 
characterized  through  normal  mode  analyses.  Energy  refinements  of  the  MP2  stractures  were  then 
performed  at  the  (JOSD(T)  level.  The  results  indicated  that  the  low-energy  path  to  decomposition  of 
syw-triazine  is  throu^  the  concerted  triple  dissociation  [Reaction  (I)],  in  agreement  with  photodissociation 
experiments  (Ondrey  and  Bersohn  1984).  Also,  the  best  prediction  of  the  decomposition  barrier 
(81.1  kcaJAnol)  is  well  below  the  experimentally  measured  upper  limit  (115  kcal/mol)  (Ondrey  and 
Bersohn  1984).  Predictions  of  the  reaction  endothermicity  of  .^/n-triazine  decomposition,  however,  were 
~8  kcalAnol  less  than  the  experimental  value  (Ondrey  and  Bersohn  1984). 

The  NDFT  results  will  be  compared  to  the  previously  reported  ab  initio  predictions  (Pai,  Chabalowski, 
and  Rice,  in  press)  and  experiment,  where  available.  In  comparing  quantum  mechanical  methods  against 
each  other  and  against  experiment,  we  will  use  the  results  calculated  with  the  cc-pVTZ  basis  set. 
Geometries,  harmonic  vibrational  fi^quencies,  and  corresponding  eigenvectors  and  relative  energies  of  the 
critical  points  will  be  compared.  Additionally,  we  will  show  the  basis  set  size  dependence  of  these 
methods. 

For  some  properties,  we  are  restricted  to  comparing  the  DFT  results  with  the  MP2  results  due  to  either 
a  lack  of  experimental  data  or  results  from  a  higher-level  correlation  technique.  In  these  comparisons,  it 
must  be  kept  in  mind  that  the  MP2  results  might  or  might  not  be  more  accurate  than  the  DFT  results. 
MP2  is  a  widely  used  method,  and  its  strengths  and  weaknesses  have  been  well-documented  for  many 
chemical  properties.  Therefore,  MP2  serves  as  a  useful  metric.  Even  with  our  rather  extensive  computer 
resources,  MP2  was  the  highest  level  of  ab  initio  theory  that  we  could  perform  on  this  size  of  system  and 
still  examine  the  basis  set  dependence. 

2.  METHODS 

DFT  calculations  of  critical  points  on  the  PES  for  jym-triazine  reactions  were  performed  using  the 
GAUSSIAN  94  program  package  (Frisch  et  al.  1995).  Geometry  optimizations  of  stable  and  transiticm 
state  species  were  done  using  the  6-31G**  (Hehre,  Ditchfield,  and  Pople  1972;  Hariharan  and  Pople  1973; 
Gordon  1980),  6-31 1++G**  (McLean  and  Chandler  1980;  Krishnan  et  al.  1980),  and  cc-pVTZ  (Woon  and 


3 


EHinning  1993;  KeiKlaU,  Dunning,  and  Harrison  1992;  Dunning  1989)  basis  sets.  The  density  functicMials 
used  are  die  various  nonlocal  spin  density  corrections  to  the  local  spin  density  (LSD)  functional  with  the 
LSD  exchange  and  the  Vosko,  Wilk,  and  Nusair  (VWN)  correlation  fimctional  (Vosko,  Wilk,  and  Nusair 
1980).  These  corrections  were  estimated  by  the  exchange  functional  of  Becke  (1988)  aiKl  the  correlation 
energy  fimction  of  Lee,  Yang,  and  Parr  (1988;  Miehlich  et  al.  1989)  (BLYP);  the  exchange  functional  of 
Becke  and  the  correlation  functional  of  Perdew  (1986)  (BP86);  the  exchange  functional  of  Becke  and  the 
correlation  functional  of  Perdew  and  Wang  (1992)  (BPW91);  and  the  exchange  functional  described  by 
the  fitted  three-parameter  hybrid  of  Becke  (1993)  and  the  correlation  functional  of  Lee,  Yang,  and  Parr 
(B3LYP).  All  geometry  optimizations  met  the  default  convergence  criteria  given  by  Gaussian  94  (Frisch 
et  al.  1995).  All  critical  points  were  characterized  through  normal  mode  an^yses.  IRC  calculations  for 
Reaction  (I)  using  all  functionals  and  the  6-31G**  basis  set  were  performed.  Ute  IRC  calculations  were 
terminated  when  minima  were  reached  as  defined  by  the  default  convergence  criteria  of  die  Gaussian  94 
set  of  programs  (Frisch  1995).  All  DFT  calculations  were  performed  using  the  default  grid  size  given  in 
Gaussian  94  (Frisch  1995). 

3.  RESULTS  AND  DISCUSSION 

3.1  Geometries.  The  geometric  parameters  predicted  at  all  levels  for  jym-triazine,  the  transition  state, 
the  weakly  bound  (HCbOj  cluster,  the  dimer  (H2C2N2),  and  the  HCN  molecule  are  listed  in  Table  1.  Ab 
initio  predictions  from  Pai,  Chabalowski,  and  Rice  (1996)  are  provided  for  comparison.  Figure  1  shows 
structures  corresponding  to  these  critical  points;  the  atom  types  and  indices  listed  in  Table  1  are  denoted 
for  jym-triazine,  (HCN)3,  HCN,  and  H2C2N2  in  Figure  1.  The  internal  coordinates  used  in  our  discussion 
are  shown  in  the  transition  state  structure  in  Figure  1(b)  or  the  (HCN)3  cluster  [Figure  1(e)].  Although 
the  atoms  are  not  labeled  on  the  transition  state  structure  in  Hgure  1(b),  the  labeling  of  the  atoms  follows 
the  same  pattern  around  the  ring  as  in  Rgure  1(a).  The  internal  coordinates  shown  in  Figure  1(b)  will  be 
discussed  later  and  are  die  same  for  all  structures  illustrated  in  Figure  1. 

3.1.1  Basis  Set  Dependence.  The  absolute  percent  differences  of  geometric  parameters  calculated 
using  the  6-31G**  or  6-31  1-h-G**  basis  sets  relative  to  those  calculated  using  die  cc-pVTZ  basis  set  are 
given  in  parentheses  within  the  columns  in  Table  1  labeled  "6-31G**"  and  "6-31 1++G**".  The  geometric 
parameters  predicted  by  all  methods  for  aU  critical  points  are  insensitive  to  increasing  basis  set  size  (<1%), 
with  two  small  exceptions.  Firsdy,  the  BLYP/6-31G**  predictions  of  the  elongated  CN  bonds  in  the 
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Table  1.  Geometric  Parameters  for  Critical  Points  on  the  5y/M-Triazine  Potential  Energy  Surface 
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Table  1.  Geometric  Parameters  for  Critical  Points  on  the  5)'m-Triazine  Potential  Energy  Surface  (continued) 


-q 

cT 


U-;  —  V*,  —  — 

-r^  —  —  —  —  —  —  —  . 


OOOOCOOOOS— “ 

rnr^tmootNao  —  O'  — 
O'/^Ot^O'Or^'O^OscT 

^0— OOOOf^f^*o 

TT-Tunoo-rroo  —  C'cMococrr 
3cr»>oo«<TooO'Ocv  —  M~,  scm 
vcro'orr'or^oor^^p.-^so 


SSS-SS-S-SSS 

'^ir^fN''j''j\OTro'0 
0'oc0^^0^0^r•,  'Tr'ir^ 
3or-icx30'0'^f-;-- 

''4vr>o<OOS?^  —  ^ 


—0-000000 


-TOt}-cxsoooO'00 

w'  *’a*  On  Os  O' 

r^w-.rsipq^g;::;;::; 


—  mOsOsoOCN'^OO 

NOsONOr^VOOOOOOOOONO-^T  — 

-TO''TO'tJ-vONONO^^^O 

—  rv,  — csj  — ppoS^S'^ 


—  -srr  —  OOOO  —  Os  —  — 
sCsnsOsoso  —  oo^x^ 
r-i  p'  oo  oo  ^  •  • 

rNjvnrJooSS);::":! 


0000000000—0 


«Nr*-rsiONr>i  —  mrsimvnv',  o 

CsOsON’^OsOOOOr'.  — 

'VO-'S-  —  '9T^r'f~'^^-_jO 

—  <<%  —  —  —  —  —  —  —  f 


rrrsj-«T  —  Tr<NTrcMr-rMCN<^ 
r-r^r-or^r^nmf'ifN4p-<N 


OOOOOOCOO 

3omoor^r-ON  —  Noirt 
aosOoooooorsir^w^JO 
r^rNt^ooeo,-:^-- 
tNvnrsippgg^^rr 


—  NO  —  OsOsOs- OO 
roosrooor~*r4ro{o 
oorMooosos^-^’- 
r^NorsiOOS'S;:;^ 


rrr^Tf  —  TfsOvONO- «ooo  — 
r-ONr^oor^TrTTTrooooNr^ 

—  cs  —  <N  —  ooogS  — 

—  fO  —  ro  —  —  —  —  —  —  —  fsi 


000—00000000 


Noosinr*-NOOs  —  onoostt  — 
Oooor-ONO'ONOP'eNHors* 
sosoNor^vor-r-'r-— -vn 
—  rsj  — cs)  — OOo2n23»n 


aorsoo  —  —  —  OvsOsO 

r^  —  r-''a’Tfo\OTr^ 

OOiOOOOsOs^-—'-' 

tN  NO  csj  o  O  S  S  ^ 


COOOOOOOO 

—  p^  —  OsOsp^moooo 

—  —  —  »r>N0O0s»O»O 

Os  NO  Os  On  On  ^  a  '  * 

r^«orJooS2);;;;r 


—  <N— <N  —  OOOO'NJ'r^NO 


r'tOscNtNr'JvOooNor^ocNso 
ON  —  qN^ososcsONTrNoooso 
'OOsoo'OP^r^p^^f^^rM 

ziSl-Si^SSSEBrS 


—0—000000 

r4  —  ps<oooO‘Nt-NOr-p«. 

sC'Ni-voooooioTrcsrM 

OsNOOsONOSrt^^"' 

raiocNppg^g-- 


NOOVOOONOP^ONp^P*.  —  —  OO 
QOONQOOOOOOOOOOOONOrsI  — 

NoosNooNor^r^r^-'^so 

—  m  —  po  —  —  —  —  —  —  —  <N 


ocMOr^p^  —  osp-r- 
so—  OOVONONOTrCsJr^ 
OOtTCOOnOn^'— -.j- 

ojNnrsjooSS^m 


o  o  o  o  o  o  o  o  o  o  oo 

wivoNONONO'OOsr'Osoooso 

—  p^  —  r"  —  osonosoocop^  — 
voONvc>o\or~p^p-,^^f^vo 

—  —  —  O,  —  qpcS^^2)q 

—  fO  —  fO  —  —  —  —  —  —  —  psj 


sssssssss 

Os^oONOpP^P~OsOs 

o  —  ooNONsocnrsJC4 

On  ■SX  On  OS  CTn  >«,•  _:  •  r 

rs^NorNipp^g-;:^ 


C'wC-w3£SSS-wC3-S 

P'  —  P-'tNP^OspMONOOOpxfO 
OsNOOs  —  Osesjcorosopsimp^ 
SONONONOVOOOOOOO^^^^ 


—0—000000 

P»*OP*OsON»o»orPTp 

lO'VNOpp  —  ooOO 

0''»'0PP^^’‘ 


NOONO  —  NOOOOsOOCOOOsOOs 
VOONNOOSNOP^P'-r'O’Nj-NONO 
NoeNNOtssovONONCi^^^rn 


o  —  o  —  OOOOOes  —  PO 


— ’pp>  —  PO  —  —  —  —  —  —  —  c 


p-ONp-NOvovoTr-'S-'ir 

tOfOPO'CTTNO'd'OOCXJ 

OOPOOOOOOO—*^-' 

<nno<npp^8-~ 


S,£,S®S,3g,o  o 

—  Os  “*  PO  fO  —  ON  O  O 
r^vop'OsONr'-rsjropn 


—  —  —  —  —  O  O  O  O  —  —  — 

cor»«ooN(OOOOOsforsirNj 
r^fop^por'ooooooos  —  ospp> 


wwSS-S-SSSS 

fOONIO»0»ONO  —  —  oo 
OsvOOsTfN3-<nvOOsoo 

aopooooooo— C<-s— I-J 
riu^pNjppg^g-^ 


_  p^  «^  — 

xuu 

2;  »- 

o  22222XXXUX2X  « 

-r*  —  pp^  Nr>  NT!  »n  p^  —  so  PM  oo  psj  E 

5  UUUUUUUU22X2  | 


—  ps*  00  Os 

U2XX 

tNrsj'(3-ooos<Nppir-k  — 

222XX2UUU 

—  pomppi  —  m-^psitN 

UUUUUU222 


Pai,  Chabalowski,  and  Rice  (1996). 
Lancaster  and  Stoicheff  (1956). 
Herzberg  (1945). 

Huber  and  Herzberg  (1979). 


reaction  [Reaction  (1)1,  (c)  the  dimer  species  H>»G<>No  associated  with  the  stepwise 
decomposition  reaction  fReaction  gpi.  (d^  the  HCN  molecule,  and  fe’t  flie  fHCNl->  cluster 
located  on  the  concerted  triple  dissociation  neacrinn  nath 


transition  state  are  2%  larger  than  the  predictions  using  the  cc-pVTZ  basis  set.  Increasing  the  basis  set 
to  6-3 1 1++G**  brings  the  BLYP  predictions  of  this  bond  lengdi  to  within  a  percent  of  the  cc-pVTZ  value. 


Secondly,  all  DFT  predictions  of  intermolecular  geometric  parameters  (N-H  distances  and  CNH  and 
NHC  angles)  in  the  weakly  bound  (HCN)3  cluster  show  a  small  basis  set  dependence  (1-5%).  Ah 
DFT/6-31G**  predictions  of  the  N-H  bond  distance  differ  by  at  least  4%  from  the  DFT/cc-pVTZ 
predictions,  but  increasing  the  basis  set  size  to  6-31 1++G**  reduces  the  difference  to  within  1%  for  all 
DFT  calculations.  Similarly,  the  DFT/6-31G**  intennolecular  bond  angles  differ  by  2-4%  from  the 
cc-pVTZ  values,  but  are  improved  upon  increasing  the  basis  set  to  6-31 1++G**.  It  is  interesting  to  note 
that  for  the  MP2  results,  increasing  die  basis  set  size  from  6-31G**  to  6-311++G**  increased  the  percent 
difference  with  the  cc-pVTZ  results  in  these  intermolecular  geometric  parameters. 
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3.1.2  Method  Comparison.  Figure  2  provides  a  comparison  of  DFT  structural  parameters  using  the 
cc-pVTZ  basis  set  with  experiment  where  available  or  with  MP2  geometries  when  experiment  is 
unavailable.  For  those  species  for  which  experimental  data  exist,  the  ordinate  axis  label  is  entitled 
"%  Error".  For  those  species  for  which  experimental  data  do  not  exist,  the  ordinate  axis  label  is  titled 
"%  Deviation"  and  refers  to  the  deviation  from  the  MP2/cc-pVTZ  predictions.  In  addition,  for  those 
structures  for  which  experimental  data  exist,  MP2  predictions  are  provided  for  additional  comparison.  A 
positive  or  negative  percent  deviation  from  the  standard  (experiment  or  MP2)  represents  an  overestimation 
or  underestimation,  respectively,  of  that  structural  parameter.  The  labels  along  the  abscissa  denote  the 
internal  coordinates  of  the  species  according  to  the  labeling  of  structures  in  Figures  1(b)  and  1(e). 


Figure  2.  Percent  error  from  experiment  for  the  structural  t?arameters  of  (a)  HCN  and  (b)  .yy/w-triazine  for 
all  methods  using  the  cc-pVTZ  basis  set  Percent  deviation  of  DFT/cc-pVTZ  methods  from 
MP2/cc-pVTZ  for  (c)  the  transition  state,  (d)  the  (HCN)3  cluster,  and  (e)  H2C2N2-  Method  is 
denoted  by  symbols  in  the  legend  in  Figure  2(a).  See  Figure  1(b)  and  (e)  for  definition  of 
internal  coordinates  shown  along  the  abscissa. 


Geometric  parameters  of  HQ^  and  aym-triazine  [Figures  2(a)  and  2(b),  respectively,]  predicted  by  all 
methods  are  within  1%  of  experiment.  DFT  predictions  of  structural  parameters  for  the  transition  state 
[Figure  2(c)]  are  within  4%  of  the  MP2  predictions,  except  for  the  CN2  bond,  which  is  being  broken  as 
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the  ^/n-triazine  proceeds  along  the  reaction  path.  For  that  parameter,  the  B3LYP  prediction  is  in  closest 
agreement  to  the  MP2  prediction.  The  BPW91  and  BP86  values  overestimate  the  MP2  value  by  ~5% 
(~0.1  A).  Similarly,  aU  of  the  DFT  predictions  of  structural  parameters  of  the  (HCN)3  cluster 
[Figure  2(d)]  are  within  3%  of  the  MP2  values,  except  for  the  weak  N-H  hydrogen  bond.  The  BP86 
prediction  of  this  bond  distance  is  in  closest  agreement  with  the  MP2  (1%),  followed  by  the  B3LYP 
prediction  (~3%).  Hie  BPW91  prediction  of  the  N-H  bond  exceeds  the  MP2  prediction  by  ~6%  (0.16  A). 
DFT  predictions  of  the  structure  of  the  dimer  species  agree  with  MP2  to  within  1.5%. 

The  good  agreement  of  DFT  with  experiment  suggests  that  the  DFT  predictions  of  stmctural 
parameters  for  stable,  covalently  btmded  structures  on  the  sym-triazine  PES  are  as  good  as  MP2.  Also, 
the  small  basis  set  dependence  of  the  DFT  geometries  indicates  that  reasonably  accurate  stmctural 
parameters  can  be  obtained  using  a  small  basis  set  (6-31G**).  The  agreement  of  DFT  with  MP2  is  better 
for  HO^,  iym-triazine,  and  the  dimer  than  for  the  transition  state  and  loosely  bound  cluster.  However, 
the  DFT  predictions  differ  from  MP2  by  no  more  than  7%. 

3.2  Vibrational  Analyses.  Normal  mode  analyses  of  aU  critical  points  on  the  sym-triazine  PES  using 
DFT  methods  and  the  6-31G**,  6-31 1++G**,  and  cc-pVTZ  basis  sets  were  performed,  and  the  harmonic 
vibrational  frequoicies  at  all  levels  are  reported  in  Table  2.  MP2  eigenvectors  corresponding  to  the 
frequencies  listed  in  Table  2  for  sym-triazine,  the  transition  state,  the  (HCN)3  cluster,  and  the  H2C2N2 
dimer  are  shown  in  Figures  3-6,  respectively. 

The  assignment  of  the  MP2  modes  for  sym-triazine  to  experiment  is  discussed  extensively  in  Pai, 
Chabalowski,  and  Rice  (1996).  The  assignments  to  experiment  are  based  upon  visual  inspection  of  the 
MP2  eigenvectors  of  sym-triazine  (Figure  3)  and  the  description  of  the  atomic  motions  given  in  the 
experimental  paper  by  Lancaster,  Stamm,  and  Colthup  (1961).  DFT  mode  assignments  are  based  upon 
comparison  of  the  DFT  eigenvectors  with  the  MP2  eigenvectors  rather  than  simply  matching  the  harmonic 
frequencies.  The  comparison  is  accomplished  by  projecting  the  DFT  eigenvectors  onto  all  choices  of  MP2 
eigenvectors  subject  to  similarity  jq)parent  from  visual  inspection  and  within  a  reasonable  energy  range. 

Assigning  some  of  the  low-fiequency  (<150  cm~^)  DFT-vibrational  modes  of  the  (HCNjj  cluster  to 
the  MP2  modes  was  difficult.  For  example.  Figure  7  shows  three  pairs  of  MP2/cc-pVTZ  and  BPW91/ 
cc-pVTZ  eigenvectors  that  we  were  unable  to  match  through  visual  inspectioa  The  MP2/cc-pVTZ 
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Table  2,  Harmonic  Vibrational  Frequencies  (cm  (continued) 
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Pai,  Chabalowski,  and  Rice  (1996). 
Herzberg  (1945). 

Hub^  and  Herzberg  (1979). 

Lancaster,  Stamm,  and  Colthup  (1961). 
Jucks  and  Miller  (1988). 


Table  2)  for  the  transition  state  for  Reaction  (I~).  Positive  and  negative  signs  denote  out-of¬ 
plane  motion. 


Table  2)  for  the  dimer.  Positive  and  negative  signs  denote  out-of-plane  motion. 
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MP2(145)  BPW91(78) 


Figure  7.  Comparison  of  MP2  eigenvectors  to  BPW91/cc-VTZ  eigenvectors  for  fa)  mode  5. 

9  of  the  (HCNV,  cluster. 


I  mode  8. 


FTiTiaCHMiiltna 


eigenvectors  correspond  to  modes  5,  8,  and  9  (See  Table  2.)  It  is  clear  that  there  are  significant 
differences  between  the  BPW91  and  MP2  eigenvectors.  The  assignments  we  made  that  are  listed  in 
Table  2  are  based  solely  on  the  combinations  that  provided  the  largest  projection  value  of  the  eigenvectors. 
For  the  MP2  and  BPW91  pairs  shown  in  Figure  7,  the  projections  of  the  BPW91  onto  the  MP2  are  0.41, 
0.74,  and  0.72,  respectively. 


In  order  to  compare  effects  of  basis  sets  and  methods  on  the  harmonic  vibrational  frequencies,  we 
calculated  the  difference  of  the  prediction  with  some  standard.  In  the  case  of  basis  set  size  dependence, 
the  standards  are  the  cc-pVTZ  calculations.  For  determination  of  the  effect  of  the  theoretical  method,  the 
standard  is  experiment,  where  available,  and  MP2  theory  where  experimental  data  do  not  exist. 
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The  theoretical  ftequencies  neglect  anhamonicity,  therefore,  the  best  theoretical  calculations  should 
in  principle  lie  above  the  exact  experimental  frequencies.  Simply  because  one  theoretical  technique  lies 
closer  to  the  exact  experimental  value  than  another  does  not  necessarily  mean  that  this  method  is  a  better 
reproduction  of  the  vibrational  modes. 

3.2.1  Basis  Set  Size  Dependence.  The  differences  in  harmonic  frequencies  calculated  using  the 
6-31G**  and  6-31 1++G**  basis  sets  relative  to  those  calculated  using  the  cc-pVTZ  basis  set  are  given  in 
parentheses  within  the  columns  in  Table  2  labeled  "6-31G**"  and  "6-31 1++G**".  DPT  predictions  of 
frequencies  of  the  critical  points  have  small  basis  set  dependencies.  Increasing  the  basis  set  from  6-31G** 
to  6-31 1++G**  brings  most  modes  to  within  ~10  cm”^  of  the  cc-pVTZ  predictions  for  all  species;  no 
frequency  differs  by  more  than  25  cm"^  In  contrast,  the  MP2  predictions  of  harmonic  vibrational 
frequencies  show  significant  basis  set  dependencies,  and  some  6-31  1-h-G**  frequencies  are  in  significantly 
poorer  agreement  with  cc-pVTZ  than  the  6-3 IG**  values. 

3.2.2  Mediod  Comparison.  Figure  8  provides  a  comparison  of  predicted  harmonic  vibrational 
frequencies  using  the  cc-pVTZ  basis  set  with  experiment  for  HCN  and  jy/n-triazine.  Figure  9  compares 
DFT-harmonic  vibrational  frequencies  with  MP2/cc-pVTZ  frequencies  for  the  transition  state,  the  (HCN)3 
cluster,  and  the  dimer. 

For  HCN  and  iyw-triazine,  BPW91,  BP86,  and  BLYP  provide  best  overall  agreement  with  experiment 
[to  within  50  cm“^  for  all  modes  except  for  sym-triazine  modes  13  [(1251  cm“^)]  and  18  [(1617  cm“^)]. 
However,  modes  13  and  18  correspond  to  two  A'2  fundamentals  which  are  inactive  and  not  observed.  The 
values  for  these  fundamentals  given  in  the  experimental  paper  are  from  "estimated  force  constants  and 
hence  are  only  very  approximate"  (Lancaster,  Stamm,  and  Coldiup  1961).  MP2  agrees  with  the 
experimental  value  given  for  5ym-triazine  mode  13  to  within  8  cm~^  However,  the  good  agreement  could 
be  fortuitous.  MP2  and  B3LYP  predictions  for  all  HCN  and  5y/n-triazine  modes  with  energies  less  than 
2000  cm“^  (except  for  mode  18,  discussed  previously)  agree  with  experiment  to  within  50  cm“^ 
agreement  with  modes  of  higher  energy  is  not  as  good. 

All  DFT-firequency  predictions  for  the  transition  state,  the  (HCN)3  cluster,  and  the  dimer  are  within 
200  cm“^  of  the  MP2  predictions.  For  bound  modes  up  to  1300  cm“^,  all  DFT  predictions  are  within 
100  cm“^  of  the  MP2  values.  B3LYP  provides  the  best  agreement  with  MP2  for  the  transition  state  and 
the  dimer.  BLYP  gives  the  best  agreement  with  MP2  for  the  (HCN)3  cluster. 
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Figure  8.  Error  from  experiment  for  the  harmonic  vibrational  frequencies  of  (a)  HCN  and  (b)  ^m-triazine 
for  all  methods  using  the  cc-pVTZ  basis  set.  Method  is  denoted  by  symbols  in  the  legend  in 
Figure  8(a). 
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Figure  9.  Difference  in  DFT/cc-pVTZ  fieouencies  from  MP2/cc-pVTZ  for  (a)  the  transition  state,  (b)  the 
(HO^^  cluster,  and  (c)  the  HoCoNo  dimer.  Method  is  denoted  by  symbols  in  the  legend  in 
Figure  9(a). 
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A  single  vibration  of  the  (HCN)3  cluster  has  been  measured  and  assigned  as  the  doubly  degenerate 
C-H  asymmetric  stretch  (Jucks  and  Miller  1988).  The  vibrational  frequency  is  3274  cm“^  and  contains 
rotationally  resolved  bands  (Jucks  and  Miller  1988).  DFT  predictions  of  the  frequency  of  this  vibrational 
mode  are  in  closer  agreement  to  erq>eriment  than  MP2.  BP86  provides  the  closest  agreement  and  is  1.3% 
larger  than  the  experiment.  Rts  of  the  bands  to  rotation^  transitions  assuming  a  planar  oblate  symmetric 
top  result  in  rotatirmal  constants  A"=B"=2C"=0.0822  cm"^  (Jucks  and  Miller  1988).  DFT  predictions  of 
the  rotaticHial  constant  deviate  from  experimait  by  1-6%,  with  BP86  and  B3LYP  having  the  best 
agreement  with  experiment  (0.0814  and  0.0815  cm~^  respectively).  MP2/cc-pVTZ  predictions  are  within 
2%  (0.0840  cm“^)  and  BLYP  and  BPW91  predictions  are  0.0791  and  0.0770  cm“^  respectively. 

The  projection  of  DFT  eigenvectors  onto  the  MP2  eigenvectors  is  not  used  only  to  assign  vibrational 
modes.  The  magnitude  of  the  projections  of  the  DFT  eigenvectors  onto  the  MP2  eigenvectors  is  a 
measure  of  the  agreement  in  the  description  of  the  atomic  motions  of  the  molecules  for  each  vibration. 
A  projection  that  has  a  magnitude  of  approximately  one  indicates  similarity  in  the  DFT  and  MP2 
predictions  of  atomic  motions  of  the  molecule  corresponding  to  that  vibrational  mode.  A  projection  that 
has  a  magnitude  near  zero  indicates  that  the  two  theoretical  methods  describe  significantly  differently  the 
motions  for  that  vibration. 

Since  the  details  of  atomic  motions  corresponding  to  each  normal  mode  of  vibration  cannot  be 
determined  through  experiment  without  assuming  some  sort  of  model  (as  in  die  experimental  work  on 
iym-triazine)  (Vosko,  Wilk,  and  Nusair  1980),  we  have  used  the  MP2  eigenvectors  as  the  standard  against 
which  to  compare  the  goodness  of  the  DFT  predictions. 

Figure  10  shows  the  projection  of  the  DFT/cc-pVTZ  eigenvectors  onto  MP2/cc-pVTZ  eigenvectors 
corresponding  to  vibrational  modes  for  sy/n-triazine,  the  transition  state,  the  (HCN)3  cluster,  and  the  dimer. 
The  eigenvectors  have  little  dependence  on  the  quality  of  the  basis  set  for  aU  of  the  critical  points  on  the 
PES.  DFT  eigenvectors  of  5ym-triazine  are  in  excellent  agreement  with  MP2  predictions  except  for 
degenerate  modes  3-4, 16-17,  and  20-21  [Figure  10(a)]. 

This  lesser  agreement  between  the  DFT  and  MP2  eigenvectors  for  these  pairs  of  degenerate  modes 
is  probably  due  to  the  nonuniqueness  of  die  eigenvectors  for  each  degenerate  pair  which  arises  from  the 
three-fold  rotational  symmetry  of  the  molecule.  Such  a  comparison  for  degenerate  modes  is  valid  only 
if  some  attempt  is  made  to  rotate  the  degenerate  DFT  eigenvectors  so  as  to  maximize  their  agreement  with 
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Figure  10.  Projection  of  DFT/cc-pVTZ  ei  genvectors  corresponding  to  harmonic  frequencies  onto  MP2/cc- 
pVTZ  eigenvectors  for  (a)  .yy/w-triazine.  (b)  the  transition  state  for  Reaction  (D.  (c)  the  (HCN)^ 
cluster,  and  (d)  the  dimer.  Symbols  in  the  legend  shown  in  Figure  10(a)  denote  the 

method. 


the  MP2  eigenvectors.  Due  to  the  good  agreement  with  the  nondegenerate  modes,  we  assumed  that  the 
lack  of  exact  agreemoit  for  these  modes  is  most  likely  an  artifact  of  the  "nonuniqueness"  of  the 
eigenvectors  for  these  double  degenerate  modes  rather  than  any  differences  in  the  predictive  quality  of 
either  theory.  Thus,  we  did  not  attempt  to  rotate  the  degenerate  DFT  eigenvectors  into  coincidence  with 
the  MP2  eigenvectors.  Excluding  the  degenerate  modes,  aU  of  the  DFT  predictions  are  in  equally  good 
agreement  with  the  MP2  vibrational  eigenvectors  of  Aym-triazine. 

Excellent  agreement  exists  among  all  DFT  and  MP2  calculations  for  modes  1-10  of  the  transition  state 
[Figure  10(b)]  as  evidenced  by  the  magnitude  of  the  DFT  projections  for  those  modes.  However,  the 
magnitudes  of  the  projections  for  the  higher  frequency  modes  of  the  transition  state  are  not  as  large,  and 
have  values  as  low  as  0.7.  fri  those  cases,  the  atomic  motions  predicted  by  DFT  for  those  vibrational 
modes  are  not  the  same  as  the  MP2  predictions,  although  the  DFT  frequencies  are  within  10%  of  the  MP2 
values.  Conversely,  the  DFT  predictions  of  the  imaginary  frequency  of  the  transition  state  ranged  from 
10  to  30%  disagreement  with  the  MP2  value  while  the  magnitude  of  the  projections  of  the  DFT 
eigenvectors  on  the  MP2  eigenvector  for  this  mode  are  all  ~1.  It  aipears  that  ^  of  the  DFT  methods 
perfomi  equally  well  in  predicting  the  vibrational  eigenvectors  for  the  transition  state. 
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For  the  (HCN)3  cluster,  the  BP86  eigenvectors  have  the  largest  projections  onto  the  MP2  eigenvectors 
[Figure  10(c)]  over  the  entire  spectrum.  Significant  scatter  exists  in  the  projectiwis  for  the  remaining 
methods,  with  BPW91  showing  the  most  notable  disagreement  with  MP2  (see  Figure  7  and  previous 
discussion).  Finally,  for  the  dimer,  all  DFT  eigenvectors  project  almost  equally  well  onto  the  MP2 
eigenvectors  over  the  entire  spectrum  [Figure  l(Kd)]. 

3.3  IRC.  Figure  11  shows  DFT/6-31G**  and  MP2/6-31G**  energies  along  the  reaction  coordinate 
for  concerted  triple  decomposition  of  5y/n-triazine  [Reaction  (I)].  Figure  12  shows  corresponding 
geometric  parameters  along  the  path.  The  aiergies  in  Figure  11  are  relative  to  three  isolated  HCN 
molecules.  Zero-point  corrections  are  not  included  in  Figure  11.  The  jyw-triazine  minimum  is  located 
on  tire  right  side  of  Figure  11,  and  the  (HCN)3  cluster  minimum  is  on  the  left-hand  side.  Although  we 
conducted  a  thorough  search,  we  were  unable  to  find  a  transition  state  connecting  minima  for  the  (HCN)3 
cluster  and  isolated  HCN  molecules. 


Figure  11.  Enereies  along  the  reaction  path  for  Reaction  <11  determined  fiom  1RC/6-31G**  calculations 
using  all  methods.  Every  tenth  point  along  the  reaction  path  from  the  IRC  calculations  is 
shown.  In  addition  to  the  (HCN)3  cluster  (illustrated  at  the  leftmost  portion  of  this  figure), 
the  transitimi  state  structure,  and  sy/n-triazine  (illustrated  at  the  rightmost  portion  of  the 
figure),  three  additional  structures  along  the  reaction  path  have  been  shown  to  enable  the 
reader  to  visualize  the  mechanism  of  concerted  dissociation  of  sy/n-triazine.  The  method  is 
denoted  by  the  symbols  defined  in  the  legend.  Energy  is  relative  to  three  isolated  HCN 
molecules. 
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Figure  12.  Internal  coordinate  changes  along  the  reaction  path  for  Reaction  (D.  See  Figure  1(b)  for  tte 
definition  of  internal  coordinates.  Bond  lengths  arc  in  A  and  bond  angles  are  in  degrces.  The 
reaction  coordinate  is  in  mass-weighted  units  of  amu^^-bohr.  Method  is  denoted  by  the 
symbols  defined  in  the  legend. 


The  energy  profiles  of  the  reaction  paths  calculated  by  all  the  methods  have  the  same  shape  but  differ 
in  the  magnitude  of  the  barrier  and  relative  well  depths  of  jym-triazine  and  the  (HCN)3  cluster. 
B3LYP/6-31G**  predicts  the  largest  classical  barrier  to  decomposition  of  sym-triazine  (94.1  kcal/mol),  and 
BLYP/6-31G**  predicts  the  smallest  (81.0  kcal/mol).  The  BP86  and  BPW91  values  (89.6  and 
90.2  kcal/mol,  respectively)  arc  within  3  kcal/mol  of  the  MP2/6-31G**  prediction  (87.2  kcal/mol).  For 
the  reverse  reaction,  the  association  of  three  HCN  molecules  to  fonn  jym-triazine,  MP2  predicts  the  largest 
classical  barrier  (47.4  kcal/mol)  and  BP86  predicts  the  smallest  (23.3  kcal/mol). 


Our  calculaticms  indicate  that  the  structures  maintain  three-fold  symmetry  throughout  the  reaction  path, 
as  can  be  seen  in  the  series  of  snapshots  of  the  atoms  at  points  along  the  reaction  path  shown  in  Figure  11. 
The  changes  in  internal  coordinates  for  all  methods,  illustrated  in  Figure  12,  arc  quantitatively  similar. 
Therefore,  DFT  predictions  of  geometric  changes  along  the  reaction  path  are  as  good  as  MP2,  and  energy 
profiles  arc  in  qualitative  agreement  with  MP2. 

3.4  Energetics.  Table  3  lists  the  DFT  and  ab  initio  absolute  and  relative  energies,  and  Table  4  lists 
DFT  and  ab  initio  energy  differences  of  critical  points  (zero  point  corrected)  on  the  potential  energy 
surface  for  sym-triazine  reactions.  The  deviation  of  eneigy  differences  calculated  using  the  6-31G**  or 
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Table  3.  Absolute  Energies  (Hartrees)  of  Species  on  the  sym-Tnazim  Potential  Energy  Surface 
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Pai,  Chabalowski,  and  Rice  (1996). 
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Pai,  Chabalowski,  and  Rice  (1996). 


6-31 1++G**  basis  sets  from  those  calculated  using  the  cc-pVTZ  basis  set  are  given  in  parentheses  within 
the  coliunns  in  Table  4  labeled  "6-3 IG**"  and  "6-311-h-G**."  Figure  13  is  an  eneigy-level  diagram 
based  on  the  values  given  in  Table  4.  This  diagram  illustrates  the  two  assumed  decomposition 
mechanisms  for  sym-triazine.  The  minimum  for  the  isolated  HCN  molecules  is  located  at  the  center  of 
the  figure  and  is  labeled  "3  HCN."  Critical  points  to  the  left  of  this  minimum  represent  features  alcmg 
the  step-wise  decomposition  pathway  [Reaction  (II)].  Critical  points  to  tte  right  of  the  minimum  labeled 
"3  HCN"  represent  features  along  the  conceited  triple-decomposition  pathway  [Reaction  (I)]. 

3.4.1  Basis  Set  Size  Eiependence.  Except  for  the  energy  difference  between  isolated  HCN  and  the 
(HCN)3  cluster,  DFT/6-31G**  values  differ  from  the  DFT/cc-pVTZ  values  by  4,4-12.4  kcalAnol.  The 
DFT/6-3 1 1-I-+G**  differences,  however,  are  within  1 .2  kcal/mol  of  the  DFT/cc-pVTZ  differences.  For  the 
reaction  barriers,  the  differences  between  DFT/6-3  11-h-G**  and  DFT/cc-pVTZ  are  no  greater  than 
0.6  kcal/mol.  While  both  the  6-311-f-K}**  and  the  cc-pVTZ  basis  sets  are  considered  triple  zeta  in  the 
valence  space,  one  very  significant  difference  between  them  is  in  the  number  of  polarization  functions  in 
each  set  For  hydrogel,  the  former  has  a  single-primitive  p-type  polarization  function,  and  a  sin^e- 
primitive  d-type  on  the  first-row  atoms.  The  cc-pVTZ  has  two  p-type  and  one  d-type  function  oa 
hydrogen,  while  the  first-row  atoms  have  two  d-type  and  one  f-type  polarization  functions,  where  aU  are 
single  uncontracted  primitive  functions.  The  DFT  results  show  very  little  dependence  on  the  quality  of 
the  basis  set  used  to  define  the  polarization  space.  Except  for  the  eneigy  difference  between  the  isolated 
HCN  molecules  and  the  (HCN)3  cluster,  ab  initio/6-3ll++G**  values  are  in  poorer  agreement  with  the 
cc-pVTZ  values  than  the  6-3 IG**  values.  For  tiie  ab  initio  calculations,  one  would  expect  the  polarization 
functions  to  play  a  significant  role  in  correcting  correlation  error  (Woon  and  Dunning  1993;  Kendall, 
Dunning,  and  Harrison  1992;  Eiunning  1989)  and  that  seems  to  be  the  case  here.  In  contrast,  the  current 
study  shows  the  DFT  results  to  have  very  little  dependence  on  the  polarization  functions. 

3.4.2  Method  Comparison.  The  zero-point  corrected  energies  of  H2C2N2  +  HCN  relative  to  three 
isolated  HCN  molecules  are  higher  than  the  zero-point  corrected  relative  energy  of  the  transition  state  for 
Reaction  (I)  for  all  methods  and  with  all  basis  sets.  Focusing  only  on  the  cc-pVTZ  basis,  the 
QCISD(T)//MP2  predicts  that  tiie  H2C2N2  +  HCN  minimum  is  6.7  kcalMol  higher  than  the  transiticxi  state 
for  concerted  triple  dissociation  [Reaction  (I)].  The  DFT  predictions  of  the  I^C2N2  +  HCN  minimum 
range  from  12.5  to  18.8  kcalAnol  higher  than  the  barrier  to  conceited  triple  dissociation.  Therefore,  evoi 
if  we  assume  that  there  are  no  barriers  to  formation  of  H2C2N2  +  HCN  from  jym-triazine  decomposition 
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Figure  13.  Biergy-level  diagram  for  the  fvm-triazine  potential  energy  surface  showing  the  minima  and 
saddle  point  energies.  Zero-point  corrected  energy  differences  between  critical  points  using 
the  different  methods  and  tte  6-31G**,  6-31 1++G**,  and  cc-pVTZ  basis  sets  are  given  in 
Table  4.  The  diagram  is  representative  of  the  features  of  the  PES  calculated  with  all  methods, 
and  is  illustrative  of  points  made  in  the  text 


or  from  association  of  two  HCN  molecules  [the  reverse  of  Reaction  (II)],  the  step-wise  decomposition 
reaction  will  be  higher  in  energy  than  the  concerted  triple  decomposition.  Our  interests  are  mainly  in  the 
low-energy  pathway  for  the  decomposition  reactions  of  sy/w-triazine;  therefore,  we  did  not  investigate  this 
step-wise  reaction  path  any  further.  Two  of  the  DFT  methods  (BPW91  and  BP86)  predict  energies  of  the 
dimer  minimum  (relative  to  isolated  HCN)  that  are  at  least  7.2  kcal/mol  lower  than  the  ab  initio 
predictions.  The  remaining  DFT  methods  are  within  2.3  kcal/mol  of  the  QCISD(T)//MP2  value. 


TTie  predictions  of  the  well  depths  of  the  (HCN)3  cluster  relative  to  the  isolated  HCN  molecules  range 
from  -4.8  to  -9.5  kcal/mol  using  the  largest  basis  set,  with  ab  initio  predicting  deeper  minima  than  DFT 
by  at  least  1.4  kcal/mol.  The  BPW91  prediction  has  the  greatest  disagreement  with  ab  initio  and  is 
4.1  kcal/mol  higher  in  energy  than  the  QCISD(T)//MP2  value.  5y/n-triazine  shows  the  opposite  trend. 
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where  BP86,  B3LYP,  and  BPW91  predict  the  relative  weU  depth  to  be  greater  than  that  predicted  by  ab 
initio  by  at  least  6,8  kcal/mol,  and  in  the  case  of  BP86  and  BPW91,  by  12-15  kcal/mol.  BLYP  prediction 
of  the  5y/n-triazine  well  depth  is  in  closest  agreement  with  the  ab  initio  predictions,  being  higher  by 
2  kcal/mol. 

With  the  exception  of  BLYP,  die  DFT/cc-pVTZ  barrier  heights  for  the  concerted  triple  decomposition 
of  jym-triazine  are  within  3  kcal/mol  of  the  ab  initio  predictions.  BPW91  and  BP86  are  within 
1.5  kcal/mol  of  the  MP2  values,  and  B3LYP  is  ~1  kcalAnol  higher  than  the  QCISD(T)//MP2  value.  BLYP 
stands  out  by  being  significantly  lower  (by  ~9-12  kcalMol)  than  jxedictions  using  all  odier  methods.  The 
only  experimental  data  for  this  barrier  were  obtained  through  fdiotodissociation  experiments,  the  results 
of  which  provided  an  ui^r  limit  of  115  kcal/mol  (Ondrey  and  Beisohn  1984). 

All  DFT/cc-pVTZ  predictions  of  the  barrier  height  for  the  reverse  reaction  of  (I),  the  concerted  triple 
association  of  3  HCN  molecules  to  form  sym-triazine,  are  lower  than  the  ab  initio  predictions  by  at  least 
4  kcal/mol.  The  B3LYP  result  is  in  closest  agreement  with  ab  initio;  all  other  methods  underestimate  the 
ab  initio  values  by  at  least  8.5  kcalAnol. 

The  measured  endothermicity  of  Reaction  (I)  is  43.2  kcalAnol  (Ondrey  and  Bersohn  1984).  Figure  14 
compares  temperature-corrected  (T  =  298  K)  reaction  enthalpies  for  .^m-triazine  decomposition  calculated 
by  all  the  methods  using  the  cc-pVTZ  basis  set  The  B3LYP  prediction  has  the  best  agreement  with 
experiment  (42,1  kcal/mol),  followed  by  BPW91  (48  kcal/mol).  BP86  predicts  a  laiger  reaction  enthalpy 
(50.4  kcalAnol),  and  BLYP  significantly  underestimates  the  experiment  (33.1  kcal/mol).  MP2  and 
subsequent  QCISD(T)  refinements  predict  reaction  enthalpies  of  35.2  and  34.4  kcalAnol,  substantially 
smaller  than  the  experiment,  and  in  almost  as  poor  agreement  as  BLYP.  This  indicates  clearly  that 
B3LYP  predictions  of  the  reaction  enthalpy  for  this  system  are  superior  to  ab  initio. 

4.  CONCLUSIONS 

In  this  study,  we  have  compared  the  perfonnance  of  NDFT  (BP86,  BLYP,  B3LYP,  and  BPW91)  using 
three  basis  sets  (6-31G**,  6-311-h-G**,  and  cc-pVTZ)  with  ab  initio  methods  [MP2  and  QCISD(T)]  and 
experiment,  where  available,  for  predicting  molecular  properties  at  critical  points  along  two  decomposition 
reaction  pathways  of  sym-triazine.  The  two  decomposition  pathways  that  are  investigated  are  a  concerted 
triple  dissociation  and  a  step-wise  decomposition  involving  an  H2C2N2  intermediate.  IMFT  and  ab  initio 
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60 


sy/77-trlazine  — >  3HCN 


Figure  14.  Temperature-corrected  fT=298  K)  reaction  enthalpies  for  Reaction  (I).  TTie  quantum 
mechanical  method  is  labeled  along  the  abscissa;  aU  calculaticms  used  the  cc-pVTZ  basis  set. 
The  experimental  value  is  given  in  Ondrey  and  Bersohn  (1984). 


geometries,  harmonic  vibrational  frequencies,  and  corresponding  eigenvectors,  and  relative  energies  of 
critical  points  were  compared.  Additionally,  energies  and  structures  were  compared  along  the  IRC  for  tl^ 
triple-concerted  decomposition  of  sym-triazine  using  the  6-31G**  basis  set.  All  methods  predict  that  the 
low-energy  pathway  to  decomposition  of  .^m-triazine  is  a  concerted  triple  dissociation,  and  that  a  weakly 
bound  cyclic  (HCN)3  cluster  exists  on  the  reaction  padi.  All  methods  predict  that  Ihree-fold  symmetry 
is  maintained  at  aU  points  along  the  reaction  coordinate.  The  energy  profile  along  the  reaction  coordinate 
is  qualitatively  similar  for  all  methods;  the  main  differences  are  in  the  sizes  of  the  reaction  barriers  and 
the  well  depths  of  sym-triazine  relative  to  three  isolated  HCN  molecules. 


We  find  that  increasing  the  basis  set  size  does  not  significantly  affect  the  predictirxi  of  structural 
parameters  for  all  methods.  Geometries  of  sym-triazine  and  HCN  predicted  by  all  methods  and  for  aU 
basis  sets  are  within  2%  of  experiment  DFT/cc-pVTZ  predictions  of  the  geometries  of  the  transition  state, 


26 


the  (HCN)3  cluster,  and  the  H2C2N2  dimer  are  in  good  agreement  with  the  MP2  predictions.  The  largest 
disagreement  between  the  DFT/cc-pVTZ  and  MP2/cc-pVTZ  structures  is  <7%. 

The  frequencies  of  the  structures  on  this  PES  show  small  basis  set  size  dependence  for  all  methods. 
The  majority  of  DFT  vibrational  frequencies  calculated  using  the  6-31 1++G**  basis  set  are  within  10  cm“^ 
of  the  cc-pVTZ  values.  BPW91,  BP86,  and  BLYP  predictions  of  the  vibrational  frequencies  of  HCN  and 
.^m-triazine  are  within  50  cm“^  of  all  experimentally  observed  modes.  MP2  and  B3LYP  predictions 
perform  as  well  as  the  other  DFT  methods  for  modes  with  energies  <2000  cm“\  All  DFT/cc-pVTZ 
predictions  of  harmcmic  vibrational  frequencies  for  the  transition  state,  the  (HCNjj  cluster,  and  the  dimer 
are  within  200  cm“^  of  MP2  predictions  for  all  modes  and  within  100  cm"^  for  modes  with  energies  less 
than  1300  cm"^  Very  little,  if  any,  experimental  data  exist  for  the  transition  state  species,  the  (HCN)3 
cluster,  and  the  dimer.  If  we  then  assume  the  MP2  results  to  be  correct,  there  spears  to  be  no  one  set 
of  fimctionals  that  stands  out  as  ctmsistently  predicting  fiequencies  in  better  agreement  with  MP2  than  the 
other  functionals.  Therefore,  we  cannot  conclude  that  one  functional  better  describes  these  species  than 
another. 

The  projection  of  the  DFT/cc-pVTZ  eigenvectors  onto  MP2/cc-pVTZ  eigenvectors  shows  that  for  most 
modes  of  .^/n-triazine  and  all  modes  of  the  dimer,  the  DFT  eigenvectors  are  quantitatively  similar  to  the 
MP2  eigenvectors.  For  rroncovalently  bound  species  [the  transition  state  and  the  (HCNjs  cluster],  the  DFT 
results  show  poorer  agreement  with  MP2  eigenvectors,  with  no  obvious  relationship  between  the  degree 
of  agreement  and  specific  regions  in  the  spectrum.  In  particular,  the  poor  agreement  between  MP2  and 
DFT  (especially  BPW91)  eigenvectors  for  a  few  low-frequency  modes  of  the  (HCN)3  cluster  and  a  few 
modes  in  the  transition  state  suggests  caution  in  theoretical  treatments  of  species  with  bonds  that  are  weak 
or  have  intermediate  bonding  character. 

Concerning  the  energetics  for  this  reaction,  there  are  some  general  trends  that  appear  in  contrasting 
the  DFT  against  the  ab  initio  results.  For  the  DFT  results,  essentially  tdl  of  the  basis  set  dependence  is 
contained  in  the  basis  functions  describing  the  valence  shell,  while  the  ab  initio  results  seem  to  depend 
equally  on  improvements  in  both  the  valence  and  polarization  bases.  These  trends  are  clearly  shown  in 
Table  4,  where  the  largest  change  in  any  aE  using  DFT  between  the  6-31 1++G**  and  cc-pVTZ  basis  sets 
is  only  1.2  kcal/mol,  while  the  change  in  going  from  the  double  zeta  6-31G**  to  triple  zeta  6-31 1++G** 
is  as  large  as  13  kcal/mol.  In  contrast,  for  the  ab  initio  calculations  in  Table  4,  the  largest  difference  in 
any  aE  going  from  6-31G**  to  6-31 1++G**  is  about  3  kcalytaol,  and  going  from  6-31 1++G**  to 
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cc-pVTZ  is  ~7  kcal/mol.  The  valence  and  polarization  functions  ^pear  to  play  an  equally  important  role 
in  the  initio  calculations,  in  contrast  to  the  DFT  where  the  polarization  functions  seem  to  play  an 
essentially  insignificant  role.  This  shows  that  the  polarization  functions  are  cmcial  in  the  initio 
calculadtms  for  correcting  the  correlation  errors  present  in  the  zeroth-order  wave  function  through  the 
perturbation  theory,  but  in  DFT  the  correlation  is  included  in  the  zeroth-order  density  which  is  apparently 
adequately  described  by  the  valence  shell  basis  set. 

A  comparison  of  the  aE’s  for  the  transition  states  and  energies  of  formation  given  in  Table  4  show 
a  wide  scatter  in  the  values  predicted  by  the  various  density  functionals.  This  is  a  noteworthy  observation, 
since  it  provides  evidence  to  die  users  of  DFT  that  transition-state  energies  and  energies  of  formation  can 
be  very  dependent  upon  the  choice  of  functionals.  Fortunately,  the  qualitative  features  of  the  PES  are  the 
same  regardless  of  the  choice  of  functionals  for  this  system.  However,  this  serves  as  a  warning  that  a 
quantitative  prediction  of  a  PES  based  upon  DFT  must  be  done  after  careful  consideration  of  the  choice 
of  exchange  and  correlation  functionals.  To  this  end,  if  one  assumes  the  QCISD(T)  results  to  be  the  most 
accurate  calculations  in  this  study,  it  would  appear  that  the  B3LYP  functionals  give  the  most  accurate  DFT 
values  for  transition-state  energies  and  energies  of  formation,  with  the  exception  of  the  zero-point 
corrected  formation  energy  of  ,^m-triazine  from  3(HCN).  The  formation  energy  ot^ained  from  BLYP, 
-30  kcal/mol,  lies  closest  to  the  QCISD(T)  value  of  -32  kcal/mol,  followed  by  the  B3LYP  value  of 
-39  kcal/mol.  This  value,  however,  is  closely  related  to  the  temperature-corrected  reaction  enthalpy  for 
jym-triazine  decomposition,  which  has  been  measured.  The  QQSD(T)-temperature-coiTected  reaction 
enthalpy  is  significantly  smaller  than  the  experimentally  measured  value  of  43.2  kcal/mol  (Ondrey  and 
Bersohn  1984),  whereas  B3LYP  is  within  1.1  kcal/mol  of  the  experiment.  With  this  exception  noted,  for 
the  set  of  functionals  used  in  this  study,  the  B3LYP  appears  to  give  the  best  overall  agreement  with  the 
QCISEKT)  results  for  aE's  corresponding  to  transition  states  and  energies  of  formation. 
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NO.  OF 

COPIES  ORGANIZATION 

2  PRINCETON  COMBUSTION 
RESEARCH  LABS  INC 
ATTN  N  A  MESSINA 
MSUMMERFELD 
PRINCETON  CORPORATE  PLAZA 
BLDGIV  SUITE  119 

11  DEERP  ARK  DRIVE 
MONMOUTH  JUNCTION  NJ  08852 

3  DIRECTOR 
SANDIA  NATL  LABS 
DIVISION  8354 
ATTN  S  JOHNSTON 
PMATTERN 

D  STEPHENSON 
LIVERMORE  CA  94550 

1  BRIGHAM  YOUNG  UNIV 

DEPT  OF  CHEMICAL  ENGINEERING 
ATTN  M  W  BECKSTEAD 
PROVO  UT  84058 

1  CALIFORNIA  INST  OF  TECH 

JET  PROPULSION  LAB 
ATTN  L  STRAND  MS  233  103 
4800  OAK  GROVE  DRIVE 
PASADENA  CA  91 109 

1  CALIFORNIA  INST  OF  TCHLGY 

ATTN  F  E  C  CULICK  MC  301  46 
204KARMANLAB 
PASADENA  CA  91 125 

1  UNIV  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTMC  LAB 
P  O  BOX  1663  MAIL  STOP  B216 
LOS  ALAMOS  NM  87545 

1  UNIV  OF  CA  BERKELEY 

CHEMISTRY  DEPARMENT 
ATTN  C  BRADLEY  MOORE 
211  LEWIS  HALL 
BERKELEY  CA  94720 

1  UNIV  OF  CA  SAN  DIEGO 

ATTN  FA  WILLIAMS 
AMESBOlO 
LA  JOLLA  CA  92093 
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NO.  OF 

COPIES  QRgANEATION 


NO.  OF 

COPIES  ORgANEATIQN 


2  UNIV  OF  CA  SANTA  BARBARA 
QUANTUM  INST 

ATTN  K  SCHOFIELD 

M  STEINBERG 

SANTA  BARBARA  CA  93106 

1  UNIV  OF  COLORADO  AT  BOULDER 
ENGINEERING  CTR 
ATTN  J  DAILY 
CAMPUS  BOX  427 
BOULDER  CO  80309-0427 

3  UNIV  OF  SOUTHERN  CA 
DEPT  OF  CHEMISTRY 
ATTNRBEAUDET 

S  BENSON 

cwrmG 

LOS  ANGELES  CA  90007 

1  CORNELL  UNIV 

DEPT  OF  CHEMISTRY 
ATTN  T  A  COOL 
BAKER  LAB 
ITHACA  NY  14853 

1  UNTV  OF  DELAWARE 

CHEMISTRY  DEPT 
ATTN  T  BRILL 
NEWARK  DE  19711 

1  UNIV  OF  FLORIDA 

DEPT  OF  CHEMISTRY 
ATTN  J  WINEFORDNER 
GAINESVILLE  FL  3261 1 

3  GEORGIA  INST  OF  TECHLGY 

SCHOOL  OF  AEROSPACE  ENGINEERING 

ATTN  E  PRICE 

WCSTRAHLE 

BTZINN 

ATLANTA  GA  30332 

1  UNIV  OF  ILLINOIS 

DEPTOFMECHENG 
ATTNHKRIER 
144MEB  1206  W  GREEN  ST 
URBANAIL  61801 


1  THE  JOHNS  HOPKINS  UNIV  CPIA 

ATTN  T  W  CHRISTIAN 
10630  LITTLE  PATUXENT  PKWY 
SUITE  202 

COLUMBIA  MD  21044-3200 

1  UNIV  OF  MICHIGAN 

GAS  DYNAMICS  LAB 
ATTNGMFAETH 
AEROSPACE  ENGNRNG  BLDG 
ANN  ARBOR  MI  48109-2140 

1  UNIV  OF  MINNESOTA 

DEPT  OF  MECHNCL  ENGNRNG 
ATTNEFIETCHER 
MINNEAPOLIS  MN  55455 

4  PENNSYLVANIA  STATE  UNTV 

DEPT  OF  MECHNCL  ENGNRNG 
ATTN  K  KUO 
MMICCI 
STHYNELL 
VYANG 

UNIVERSITY  PARK  PA  16802 

2  PRINCETON  UNIV 
FORRESTAL  CAMPUS  LIB 
ATTNKBREZINSKY 
IGLASSMAN 

PO  BOX  710 
PRINCETON  NJ  08540 

1  PURDUE  UNIV 

SCHL  OF  AERO  &  ASTRO 
ATTN  JR  OSBORN 
GRISSOM  HALL 
WEST  LAFAYETTE  IN  47906 

1  PURDUE  UNIV 
DEPT  OF  CHEMISTRY 
ATTN  E  GRANT 

WEST  LAFAYETTE  IN  47906 

2  PURDUE  UNIV 

SCHL  OF  MECHANICAL  ENGNRNG 
ATTN  N  M  LAURENDEAU 
SNBMURTHY 
TSPC  CHAFFEE  HALL 
WEST  LAFAYETTE  IN  47906 
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NO.  OF 

COPIES  ORGANIZATION 

1  RENSSELAER  POLYTCHNC  INST 

DEPT  OF  CHEMICAL  ENGNRG 
ATTNAFONTUN 
TROY  NY  12181 

1  STANFORD  UNIV 

DEPT  OF  MECHNCL  ENGNRNG 
ATTN  R  HANSON 
STANFORD  CA  94305 

1  UNIVOFTEXAS 

DEPT  OF  CHEMISTRY 
ATTN  W  GARDINER 
AUSTIN  TX  78712 

1  VA  POLYTECH  INST 

AND  STATE  UNTV 
ATTNJASCHETZ 
BLACKSBURG  VA  24061 

1  APPLIED  COMBSTN  TECHLGY 
ATTN  AM  VARNEY 

PO  BOX  607885 
ORLANDO  FL  32860 

2  APPLIED  MECHANICS  REVIEWS 
ASME 

ATTN  R  WHITE  &  A  WENZEL 
345  E  47TH  ST 
NEW  YORK  NY  10017 

1  TEXTRON  DEFENSE  SYSTEMS 

ATTN  A  PATRICK 
2385  REVERE  BEACH  PKWY 
EVERETT  MA  02149-5900 

1  BATTELLE 

TWSUAC 
505  KING  AVE 
COLUMBUS  OH  43201-2693 

1  COHEN  PROFESSIONAL  SVCS 

ATTN  NS  COHEN 
141  CHANNING  ST 
REDLANDS  CA  92373 

1  EXXON  RESEARCH  &  ENG  CO 

ATTN  A  DEAN 
ROUTE  22E 
ANNANDALENJ  08801 


NO.  OF 

COPIES  ORGANIZATION 

1  GENERAL  APPLIED  SCIENCE 

LABS  INC 
77  RAYNOR  AVE 
RONKONKAMA  NY  11779-6649 

1  GENERAL  ELECTRIC  ORDNANCE 

SYSTEMS 
ATTNJMANDZY 
100  PLASTICS  AVE 
PITTSFIELD  MA  01203 

1  GENERAL  MOTORS  RSCH  LABS 
PHYSICAL  CHEMISTRY  DEPT 
ATTNTSLOANE 

WARREN  MI  48090-9055 

2  HERCULES  INC 
ATTNWBWALKUP 
E  A  YOUNT 

PO  BOX  210 
ROCKET  CTR  WV  26726 

1  HERCULES  INC 

ATTN  R  V  CARTWRIGHT 
100  HOWARD  BLVD 
KENVILNJ  07847 

1  ALLIANT  TECHSYSTEMS  INC 

ATTN  RE  TOMPKINS 
MN  11  2720 

600  SECOND  ST  NORTH 
HOPKINS  MN  55343 

1  IBM  CORPORATION 

RESEARCH  DIVISION 
ATTNACTAM 
5600  COTTLE  RD 
SAN  JOSE  CA  95193 

1  nr  RESEARCH  INST 

ATTNRFREMALY 
10WEST35THST 
CHICAGO  IL  60616 

1  LOCKHEED  MIS  &  SPACE  CO 

ATTN  GEORGE  LO 
3251  HANOVER  ST 
DEPT  52  35B204  2 
PALO  ALTO  CA  94304 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  OLIN  ORDNANCE 

ATTN  V  MCDONALD  LIB 
PO  BOX  222 

ST  MARKS  FL  32355-0222 

1  PAUL  GOUGH  ASSOCIATES  INC 

ATTN  PS  GOUGH 
1048  SOUTH  ST 
PORTSMOUTH  NH  03801-5423 

1  HUGHES  AIRCRAFT  CO 

ATTNTEWARD 
8433FALLBROOKAVE 
CANOGA  PARK  CA  91303 

1  SCIENCE  APPLICATIONS  INC 

ATTNRBEDELMAN 
23146  CUMORAH  CREST 
WOODLAND  HILLS  CA  91364 

3  SRI  INTERNATIONAL 

ATTN  G  SMITH 
D  CROSLEY 
D GOLDEN 

333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

1  STEVENS  INST  OF  TECH 

DAVIDSON  LAB 
ATTNRMCALEVYin 
HOBOKEN  NJ  07030 

1  SVERDRUP  TECHLGY  INC 

LERC  GROUP 
ATTNR  J  LOCKE  MS  SVR  2 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

1  SVERDRUP  TECHLGY  INC 

ATTNJDEUR 

2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

3  THIOKOL  CORPORATION 

ELKTON  DIVISION 
ATTNR  BIDDLE 
RWILLER 
TECH  LIB 
PO  BOX  241 
ELKTON  MD  21921 


3  THIOKOL  CORPORATION 

WASATCH  DIVISION 
ATTN  SJ  BENNETT 
PO  BOX  524 

BRIGHAM  CITY  UT  84302 

1  UNITED  TECHNOLOGIES 

RESEARCH  CTR 
ATTNACECKBRETH 
EAST  HARTFORD  CT  06108 

1  UNITED  TECHNOLOGIES  CORP 

CHEMICAL  SYSTEMS  DTV 
ATTN  RR  MILLER 
PO  BOX  49028 
SAN  JOSE  CA  95161-9028 

1  UNIVERSAL  PROPULSION  CO 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVE 
PHOENIX  AZ  85027-7837 

1  VERTTAY  TECHLGY  INC 

ATTNEBnSHER 
4845  ME.LERSPORT  HWY 
EAST  AMHERST  NY  14051-0305 

1  FREEDMAN  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  RD 
BALTIMORE  MD  21209-1525 

6  ALLIANT  TECHSYSTEMS  INC 

ATTN  C  CANDLAND 
L  OSGOOD 
R BECKER 
JBODE 
RBURETTA 
M  SWENSON 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  DIRECTOR 

US  ARMY  BENET  LABS 
ATTN  AMSTA  AR  CCB  T 
SAM  SOPOK 
WATERVLIET  NY  12189 
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NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PRQVINO  OROUNP 

35  DIRUSARL 

ATTN:  AMSRL-WT-P,  A  HORST 

AMSRL-WT-PC, 

BEFORCH 

GF  ADAMS 

WR  ANDERSON 

RABEYER 

SWBUNTE 

CFCHABALOWSKI 

K  P  MCNEILL-BOONSTOPPEL 

A COHEN 

RCUMPTON 

R  DANIEL 

DDEVYNCK 

NFFELL 

JMHEIMERL 

B  HOMAN 

AJKOTLAR 

WFMCBRATNEY 

KLMCNESBY 

NEMEGHER 

MS  MILLER 

AWMIZIOLEK 

JB  MORRIS 

JE  NEWBERRY 

S  ANEWTON 

R  A  PESCE-RODRIGUEZ 

BMRICE 

LSEGER 

PVSHARMILA 

R  C  SAUSA 

MASCHROEDER 

W  TEAGUE 

JAVANDERHOFF 

DVENIZELOS 

AWHREN 

C  WILLIAMSON 


37 


Intentionally  left  blank. 


38 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratoiy  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers  to 
the  items/questions  below  will  aid  us  in  our  efforts. 

1 .  ARL  Report  Number/Author  ARL-TR-1217  rPai') _ Date  of  Report  October  1096 _ 

2.  Date  Report  Received  _ _ _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report  will 

be  used.) _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  Gaieral  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to  organization, 
technical  content,  format,  etc.) _ _ _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFHCtAL  BUSINESS 

BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001  ,APG,MD 


POSTAGEWILL  BE  PAID  BY  ADDRESSEE 
DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRLWTPC 

ABERDEEN  PROVING  GROUND  MD  21005-5066 


